Entanglement in a Solid State Spin Ensemble 
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Entanglement is the quintessential quantum phenomenon and a necessary ingredient in most 
emerging quantum technologies, including quantum repeaters [T], quantum information processing 
(QIP) [2] and the strongest forms of quantum cryptography [3 . Spin ensembles, such as those 
in liquid state nuclear magnetic resonance [4] |5] , have been powerful in the development of quan- 
tum control methods, however, these demonstrations contained no entanglement and ultimately 
constitute classical simulations of quantum algorithms. Here we report the on-demand genera- 
tion of entanglement between an ensemble of electron and nuclear spins in isotopically engineered 
phosphorus-doped silicon. We combined high field/low temperature electron spin resonance (3.4 T, 
2.9 K) with hyperpolarisation of the 31 P nuclear spin to obtain an initial state of sufficient purity 
to create a non-classical, inseparable state. The state was verified using density matrix tomography 
based on geometric phase gates, and had a fidelity of 98% compared with the ideal state at this field 
and temperature. The entanglement operation was performed simultaneously, with high fidelity, 
to 10 10 spin pairs, and represents an essential requirement of a silicon-based quantum information 
processor. 



Most QIP algorithms applied to spin ensembles have 
been implemented in a regime of weak spin polarisation. 
However, due to the very low purity of the states used, 
any exponential enhancement offered by quantum me- 
chanics disappears when the scaling of total resources 
is considered. Highly mixed, or weakly initialised, en- 
sembles are often interpreted as the sum of a perfectly 
mixed component (given by a normalised identity ma- 
trix in the density matrix representation) and a smaller 
pure component: p true = (1 — e)I/d+ epo, where d is the 
dimensionality of the state. The I component is invari- 
ant under unitary operations and not directly observable 
by magnetic resonance, which produces measurements of 
the population differences across allowed electron- and 
nuclear-spin transitions. It is therefore straightforward 
to ignore the maximally mixed component: this approach 
is called the "pseudo-pure approximation" [6]. 

There are a number of entanglement witnesses or 
monotones which are able to distinguish entangled states 
from (classical) separable ones [7]. A widely used test is 
the positive partial transpose (PPT) criterion, which is 
both a necessary and sufficient test of entanglement for 
two coupled spin- 1/2 particles [EJ[9]. Applying this test 
to the mixed state p t rue above, it can be easily shown 
that the minimum value of e which permits the overall 
state to be entangled is 1/3 [4]. 

Typical values for e in liquid state NMR and ESR (at 
X-band, 10 GHz, 5 K) are e ~ 10" 5 and e ~ 10~ 2 , respec- 
tively, well below the required threshold for the PPT test. 
Thus, while experiments performed in this regime pro- 
vide a valuable test-bed for techniques in entanglement 



generation and detection [5j[T2j[T3], the states created are 
only pseudo-entangled, and fully separable. (A notable 
exception was the use of chemical methods to generate 
highly polarised hydrogen spin pairs [14] , though this is a 
single-shot experiment with limited scalability.) In order 
to overcome this limit, we require states of higher initial 
purity, and a method to measure the I component of the 
density matrix. 

We follow a hybrid approach, using both the electron 
and nuclear spin associated with a phosphorous donor in 
silicon. Isolated donors in isotopically engineered semi- 
conductors are of particular interest as they possess ex- 
cellent decoherence characteristics (both electron and nu- 
clear T2 exceed seconds [I5j [16] ) , can be controlled with 
high-fidelity using microwave and radio-frequency pulses 
[T71 [T8] and are promising for integrating quantum tech- 
nologies into conventional semiconductor devices [19]. 

Neglecting the weak polarisation of the nuclear spin, 
the initial state populations are determined by the elec- 
tron spin Zeeman energy, as shown in Figure [TJa) where 
a = exp(— q\±bB jkBT\ g is the electron g-factor, \±b is 
the Bohr magneton, /c# is Boltzman's constant and B 
and T are the experimental magnetic field and temper- 
ature respectively. At a high magnetic field (3.4 T) and 
low temperature (2.9 K), the donor electron spin is ther- 
mally polarised to ~ 66%, however the 31 P nuclear spin, 
with a much weaker magnetic moment, has only ~ 0.04% 
polarisation. Various methods exist for indirectly trans- 
ferring electron spin polarisation to the nuclear spin, un- 
der the heading dynamic nuclear polarisation [20j [21], 
and often exploiting cross-relaxation process involving 
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FIG. 1. Sequences for nuclear spin hyperpolarisation and 
entanglement generation for this coupled S = 1/2, / = 1/2 
spin system, a) The initial state is at thermal equilibrium 
where populations (shown in green) are distributed accord- 
ing to the electron spin polarisation at this magnetic field 
and temperature (see text). A pair of microwave and ra- 
diofrequency tt pulses (tt 1,3 and tt 3,4 ) are applied to move 
spin populations to favour the |f) nuclear spin state. After 
some time twait, much longer than the electron spin relax- 
ation time (Tie), there is significant majority population in 
the |ti) state. Nuclear spin and cross-relaxation processes 
are much slower than Ti e . b) Illustration of the 28 Si:P single 
crystal host with a substitutional phosphorus donor and its 
electron spin, c) Starting from the hyperpolarised state in 
(a), an electron spin coherence is generated by a microwave 
(tt/2) 1,3 pulse. A radiofrequency tt 3,4 pulse transforms this 
into the final entangled state, containing a superposition of 
|tt) an d \H}- d) The growth in the electron spin echo inten- 
sity measured on the |1) : |3) transition provides a measure 
of the population ratio a. e) This hyperpolarisation sequence 
minimises the linear entropy of the two-spin state for a given 
value of a. 



simultaneous electron and nuclear spin flips. Here, we 
exploit the relative absence of cross-relaxation leading 
to a substantial difference in the relaxation times of the 
electron and nuclear spin [17 in order to hyperpolarise 
the nuclear spin rapidly and with high efficiency. This 
hyperpolarisation process is similar to "algorithmic cool- 
ing" methods, where a particular qubit relaxes quickly 
due to coupling to a heat bath [3]. 

Figure [TJa) illustrates our method for tackling the twin 
challenges of measuring and minimising the I compo- 
nent in the density matrix of the coupled electron- nuclear 
spin system. The hyperpolarisation of the nuclear spin 
can be understood as a SWAP operation with the (ther- 
mally polarised) electron spin, using a combination of 
resonant microwave (MW) and radiofrequency (RF) tt 
pulses. This is followed by a delay twait which is sub- 
stantially longer than the electron spin relaxation time 
Ti e (specifically twait ~ 8Ti e ), during which the elec- 
tron spin relaxes back to thermal equilibrium. On this 
timescale, other relaxation processes (such as pure nu- 



clear spin flips, or electron- nuclear spin flip-flops) are 
orders of magnitude slower and can be neglected. The 
resulting hyperpolarised state is: 

p = ^(a\l)(l\ + a 2 |2)(2| + |3> <3| + a|4)(4|), (1) 

where Z = 2(1 + a) is a normalising constant. 

Although spin-echo sequences can only be used to 
probe the population differences across energy levels, we 
can obtain a direct measure of the population ratio a 
by measuring the electron spin-echo amplitude between 
levels |1) and |3) before and after the hyperpolarisation 
sequence, as shown in Figure [TJd). Due to the enhanced 
polarisation of the nuclear spin, a spin-echo measured on 
this transition increases by a factor 2/(a + 1), compared 
to the measurement from a fully-relaxed thermal state. 
This measure is strictly conservative: it places a lower 
bound on the true polarisation of the electron as imper- 
fections such as pulse errors or residual relaxation pro- 
cesses only lead to a lower apparent state purity. Using 
the measure, we observe an enhancement of the echo in- 
tensity by a factor of 1.643(2), corresponding to a upper 
bound of a < 0.217(2). 

This hyperpolarisation sequence corresponds to a de- 
crease in linear spin entropy, made possible by the open 
quantum system's contact with the lattice heat bath (see 
Figure [lje)). Importantly, this approach leads to the 
minimum possible linear entropy given the electron spin 
polarisation resource and type of relaxation present [3]. 
Entanglement is maximised in a mixed 2-qubit density 
matrix by first minimising the linear entropy, and then 
generating an entangled coherence across the levels with 
the greatest and second- smallest population [4j [23] . Fol- 
lowing this strategy, we create an entangled state using a 
coherence-generating MW ^ 1,s pulse, followed by an RF 
7r 3 ' 4 pulse (as shown in Figure [ljb)), yielding the target 
state: 
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This density matrix is entangled according to the PPT 
criterion when a < 0.432, while other preparation meth- 
ods (such as pseudo-pure state preparation) require sub- 
stantially higher polarisation (see Supplementary Infor- 
mation) . 

Having prepared the initial state and performed an en- 
tangling operation, we now use density matrix tomogra- 
phy to extract the final two-spin state. Due to the weak 
magnetic moment of nuclear spins and necessarily low 
donor concentration in our sample, we are restricted to 
non-projective measurements of the electron spin ensem- 
ble along the a x and a y bases which can be performed 
selectively on the mi state of the nuclear spin (in prod- 
uct operator formalism, these bases can be written as 
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FIG. 2. Electron/nuclear spin phase rotations reveal the off- 
diagonal elements of the density matrix, a) Under the appli- 
cation of two consecutive tt 1,3 pulses around different axes (0), 
the eigenstates |1) and |3) undergo closed trajectories on the 
Bloch sphere with equal and opposite solid angle £1 — ±20. 
Each state picks up a phase equal to half this solid angle, b) 
This 7To — phase gate is applied to both electron |1):|3) and 
nuclear |3):|4) transitions, where the two phases are varied by 
different increments 5<j) and 5 a as the experiment is repeated. 
Example oscillations are shown for three experiments where 
we generate an electron coherence |1)(3|, nuclear coherence 
|4)(3| and double quantum coherence |2)(3|. c) Fourier trans- 
forms of the oscillations with respect to increment number 
show peaks located at 0.050(8), 0.031(5) and -0.079(8), in ex- 
cellent agreement with the set frequencies v<$> — 2n /5<j) = 0.05 
and v a 2tx/5cj 0.03. 



Diagonal elements of the density matrix (correspond- 
ing to state populations) are obtained by mapping pairs 
of population differences into an electron spin echo on 
the |1):|3) transition (S Xjy I a ). The accurate detection 
of off-diagonal elements (coherences) is a more elabo- 
rate process, made by selectively labelling the coherence 
between each pair of eigenstates with a distinguishable 
time- varying phase [5[|24]. Under this process, a partic- 
ular phase-accumulation rate provides the signature of a 
particular coherence, allowing the off-diagonal elements 
to be reconstructed from the amplitudes in the Fourier 
transform of a measured signal. 

Here we follow an approach inspired by the Aharonov- 
Anandan geometric phase gate [6j [7] to apply arbitrary 
phases in a fixed time to the four different eigenstates, 
and thus separately label each of the possible coherences. 
We apply two tt pulses, along different axes, across a 
transition between a pair of eigenstates. The phase ac- 
quired by each eigenstate is opposite and equal to half 
the solid angle of its trajectory on the Bloch sphere (see 
Figure 2[ a)). Thus, applying 7Tq' 3 followed by — 7r^' 3 (sub- 
scripts aenote pulse phase and thus nominal rotation 
axis) leads eigenstates |1) and |3) to undergo trajecto- 
ries of equal and opposite solid angle ±20. A similar 
operation is applied to the nuclear spin transition: 7Tq' 4 , 
— 7r 3 ' 4 , such that the total operator describing the action 



element being read 

echo 




FIG. 3. Measuring an entangled density matrix, a) The 
full pulse sequence used to prepare, entangle and measure the 
two-spin state. The final readout stage was changed according 
to the density matrix element being measured: two examples 
are shown for the |1)(2| and |1)(4| states, b) The obtained 
density matrix is shown as solid bars, while the dashed outline 
(zero where not shown) shows that of an ideal state given 
a = 0.217. The fidelity of the ideal state with the measured 
density matrix is 98%. 
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The value of <j) is incremented by 5(j) on each shot of 
the experiment, with effective frequency v$ = 2ir/5(j) 
(and similarly for a, 5 a and u a ). We then map each 
off-diagonal element of the density matrix in turn into 
S Xj yI a using a set of appropriate MW and RF tt pulses, 
and measure the amplitude of the Fourier component at 
the effective frequency corresponding to that coherence. 
Quadrature measurement allows us to discriminate be- 
tween positive and negative frequencies. The presence of 
other Fourier peaks would be illustrative of pulse errors in 
the mapping sequence, but as seen in Figure [2|c), such 
errors are negligible, even in the absence of operations 
such as phase cycling. 

Combining our measurements of the identity compo- 
nent, diagonal, and off-diagonal elements of the density 
matrix of the electron-nuclear spin system, we obtain: 
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This state has a minimum eigenvalue under the PPT 
test of —0.19(1) and a concurrence (C) of 0.43(4), each of 
which confirm the presence of finite entanglement. The 
fidelity of the measured density matrix with the target 
state given a = 0.217 is 98.2(2)%, and results of this to- 
mography process are shown in Figure [3j To obtain the 
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uncertainty in these values we used Monte Carlo gener- 
ation of physical density matrices based on the standard 
error of each matrix element due to noise (see Supple- 
mentary Information). 

The finite entanglement shown can offer direct advan- 
tages over classical methods in applications such as quan- 
tum sensors [271 [28] . To achieve higher purity entan- 
gled states one could move to lower temperatures e.g. we 
would expect C ~ 0.99 if these experiments were per- 
formed at 0.8 K. Complementary to this approach, en- 
tanglement purification could be performed using a larger 
Hilbert space at each node [29] , for example using a donor 
atom with a higher nuclear spin (such as bismuth with 
7 = 9/2 [301131]). 

The electron-nuclear spin entanglement generated here 
could also be mapped into an entangled state between nu- 
clear spin pairs. By SWAPping the state of the electron 
spin with a second, coupled nucleus, for example, one 
could attain nuclear-spin entanglement in a regime where 
the thermal polarisation of the nuclei would be orders 
of magnitude too small and the direct coupling between 
them weak. Clusters of up to eight nuclei coupled to a 
single electron spin have been explored in other materi- 
als [32], though the scaling of such an approach appears 
limited. A scalable network of entangled nuclear spins 
could be generated by exploiting the ability to ionise 
the donor and transfer the electron onto a neighbour- 
ing donor site [33j [34] . These operations, combined with 
single-shot readout of the P-donor spin [35] and globally- 
controlled electron- nuclear spin entanglement such as we 
have demonstrated, form the basis for a cluster-state 
quantum computer based in silicon [36] . 
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Supplementary information 
Materials and Methods 

Si:P consists of an electron spin S=l/2 (g = 1.9987) coupled to the nuclear spin / = 1/2 of 31 P through an isotropic 
hyperfine coupling of a = 4.19 mT [1]. The W-band EPR signal comprises of two lines (one for each nuclear spin 
projection Mj = ±1/2). Our experiments were performed on the low-field line of the EPR doublet corresponding to 
Mi = 1/2. At 2.9 K and 3.36 T, the electron and nuclear spin T\ were measured to be approximately 0.6 s and 100 s, 
respectively. 

The sample consists of a 28 Si-enriched single crystal about 0.5 mm in diameter with a residual 29 Si concentration 
of order 70 ppm, produced by decomposing isotopically enriched silane in a recirculating reactor to produce poly-Si 
rods, followed by floating zone crystallisation. Phosphorus doping of ~ 10 14 cm -3 was achieved by adding dilute PH3 
gas to the Ar ambient during the final float zone single crystal growth. Further information on the sample growth 
has been reported elsewhere [2]. 

Pulsed EPR experiments were performed using a W-band (94 GHz) Bruker Elexsys 680 spectrometer, modified to 
allow microwave phase control, and equipped with a 6 T superconducting magnet and a low temperature helium-flow 
cryostat (Oxford CF935). The cryostat was pumped to achieve a temperature of 2.88 K (internal thermocouple) with 
a consistent upper temperature limit confirmed by the spin temperature measurement (see text). Typical pulse times 
were 56 ns for a mw tt pulse and 100 /is for an rf tt pulse. To achieve arbitrary phase control, RF pulses were generated 
using a Rohde and Schwarz AFQ100B together with an Amplifier Research 500 W amplifier. 



Hyperpolarisation 

Si:P is described by an isotropic spin Hamiltonian (in angular frequency units): 

Ho = oj e S z - ujiI z + a-S-I, (4) 

where u e = gf3Bo/h and ojj = gif3 n B /h are the electron and nuclear Zeeman frequencies, g and gi are the electron 
and nuclear g-factors, j3 and f3 n are the Bohr and nuclear magnetons, H is Planck's constant and Bo is the magnetic 
field applied along z-axis in the laboratory frame. S XjVjZ = <J x ,y,z ® I2 and I x ,y,z = I2 ® &x,y,z- 

In the spin basis introduced in the main text, the implementation of the hyperpolarisation sequence upon the 
first-order Boltzmann thermal state is given by: 



e — huj e S z /k B T 



(5) 



= |(a|l)(l|+a|2)(2| + |3)(3| + |4)(4|) 
4' 3 -> P =|(|1)(1| + «|2><2| + a|3>(3| + |4)(4|) 
TTo' 4 P =^(|1)(1| + a|2)<2| + |3)(3| + a|4)(4|) 
twait -»• P ~{a\l){\\ + a 2 |2)(2| + |3> <3| + a|4)<4|) 

where a = e - huJ e /k B T an( ^ g = 2(1 + a), T is temperature in Kelvin and co e the resonant frequency of the electron 
spin. The ratio of spin-echo amplitudes between the hyperpolarised state to the thermal state on the |1):|3) transition 
is therefore given by 4/Z — 2/(a + 1). 

This hyperpolarisation process minimises the linear entropy of the system subject to the available physical pro- 
cesses [3]. Many measures can be used to characterise a state's mixedness; in this work we will use linear entropy 
because its upper limit for entangled states is tighter than the von Neumann entropy limit [4J. The linear entropy is 
given by 

^W^^-Il- traced 2 )] (6) 
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and so the linear entropies of the thermal and hyperpolarised states are given by 

2(1 + 4a + a 2 ) 



SlXPthermal) — 
Sl (Phyperpol) 



3(1 + a) 2 
16a(l + a + a 2 ) 



3(1 + a) 4 
The linear entropy decrease is: 



<7 (n \ q (n \ 2(l-a) 2 (l + a 2 ) 

^MPthermalJ ~ (Phyperpol ) = 3^ + a y V'J 



Entanglement Thresholds 

As discussed in the text, our target entangled state was optimally chosen within the limits of the state's linear 
entropy. We can apply the entangling sequence to the hyperpolarised state expressed in the form of a density matrix 
in the {|1),|2),|3),|4)} basis. 
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1 -a\ 
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2a 2 


















2a 







I 1 


— a 








1 + a/ 



(8) 



The partial transpose of this matrix has a zero crossing when 

a 3 -(l-a) 2 /4 = 

or when a w 0.432. 

It is worth briefly considering the crossing-points for other preparation strategies, in each case using the optimal 
entangling sequence described above. Three examples include i) the pseudopure preparation scheme applied to the 
thermal state ii) the thermal state with no preparation and iii) the pseudopure preparation scheme applied to the 
hyperpolarised state. For i), the prepared pseudopure state has the form 



1 



/(l + 2a)/3 















(l + 2a)/3 
1 
(l + 2a)/3/ 



This initial state produces the maximal amount of entanglement when arranged in the form equivalent under local 
unitaries to 



1 



/(4 + 2a)/6 

(1 + 2a)/3 

\(2 - 2a)/6 



(2 - 2a)/6\ 

(l + 2a)/3 

(4 + 2a)/6/ 



The partial transpose of this matrix has a zero crossing when 



2a + a 2 







and so the only 'physical' zero crossing occurs at a = 0. Of course, these calculations assume that the nuclear spin 
polarisation is negligible; in the regime of a « the polarisation of the nucleus would have to be considered in this 
calculation. This would allow for a small but nonzero threshold for a, dependent upon the nuclear spin isotope. 
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For ii), the maximally entangled state is equivalent under local unitaries to 



Z 



/(l + a)/2 (l-a)/2\ 

a 

1 

\(l-a)/2 (l + a)/2/ 



The partial transpose of this matrix has a positive zero crossing when a = — 3 + 2V% so entanglement is possible 
below a « 0.17. 

For iii), the maximally entangled state is equivalent under local unitaries to 



'Z 2 



/(3 + 2a + a 2 )/6 (3 - 2a - a 2 )/6\ 

a(2 + a)/3 

a(2 + a)/3 

\(3-2a-a 2 )/6 (3 + 2a + a 2 )/6/ 



The partial transpose of this matrix has a positive zero crossing when a = y/2 — 1, and so entanglement is possible 
below a « 0.4142. 



Density Matrix Calculations 

We used phase (Z) rotations to probe the existence of coherences without disturbing the state populations, which 
were measured separately by mapping pairs of population differences into the S x observable. In typical electron spin 
resonance experiments, phase rotations are not directly available operations; they can be implemented with various 
strategies including the use of off-resonant pulses and appropriate delays, phase-shifted final reference frames [5] and 
composite rotations made up of pure X and Y rotations of variable length [6]. In this work, we follow an approach 
inspired by the Aharonov-Anandan geometric phase gate [7] and apply two selective MW or RF tt pulses of differing 
phases. Given the larger Hilbert space of this two-spin system, and the fact that we are applying selective rotations 
that cannot be understood in terms of the independent manipulation of either spin, the most instructive way to 
appreciate the effect of this gate is to examine the phase acquired by each eigenstate under the rotation (see Figure 
2 of the main text). However, for a single spin 1/2, the same gate has a convenient visual representation that is 
illustrated in Figure [4] in terms of the rotation of the Bloch sphere. 

To calculate the final density matrix we first constructed the pseudo-pure density matrix. A baseline measurement 
was taken as an average of 2000 samples, and all datasets were baseline-corrected before processing. The population 
differences were measured by an average of 100 samples and scaled with respect to a measured thermal amplitude 
(also taken as an average over 100 samples), and adjusted to have unit trace with the addition of an appropriately 
scaled identity matrix. 

The coherence-measurement pulse sequences were all individually tested by verifying their output frequencies when 
applied to their target coherences. Upper diagonal elements were gathered and their conjugate values populated the 
lower diagonal elements. The coherences were collected with 128 points and baseline-corrected before being Fourier 
transformed and normalised with respect to the thermal amplitude. A narrow integral over the appropriate frequency 
position was made to measure the coherences. The frequencies were chosen such that all frequency peaks were 
well-resolved with a 128-point Fourier transform. 

The calculated pseudopure matrix p pp was added to the appropriate amount of identity matrix I as determined by 
a spin-temperature measurement, where e - huj s/k b T = a < 0.217. The explicit reconstruction is given by 

p F = [a/ (2(1 + a))]I + [(1 - a)/((l + a))]p pp 

The fidelity of the final density matrix pp with pp, the target density matrix (calculated by Equation [8] with 
a = 0.217), was calculated as 




as proposed elsewhere [8]. 

The errors corresponding to each element were calculated according to the standard error of the measurement taken. 
For the populations this consisted of the standard error of the direct difference measurements; for the coherences this 
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# 

Z(0) Z{n/2) Z(?r) 

Effective Rotation 

FIG. 4. Geometric Phase Rotations Geometric Z-gates can be applied to a single spin using two consecutive tt pulses along 
different axes. The black arrows at each step illustrate the applied rotation to the Bloch sphere, and the colouring displays this 
rotation. Three sample geometric Z-gate rotations are shown along the bottom. 

consisted of the standard error of the Fourier-transformed signal over the integral peak width. These density matrix 
element errors were transformed into final negativity, concurrence and fidelity errors by Monte Carlo generation of 
density matrices. The generated matrices deviated from the measured matrix in each element by an amount chosen 
randomly from a normal distribution whose standard deviation matched that elements' error. Once re- normalised, 
unphysical matrices were discarded and statistics on physical matrices were collected. In total, 2 12 matrices were used 
to compile the final errors. 
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